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The close theoretical analogy between the physics of rapidly rotating atomic Bose condensates 
and the quantum Hall effect (i.e. a two dimensional electron gas in a strong magnetic field) was first 
pointed out ten years ago. As a consequence of this analogy, a large number of strongly correlated 
quantum Hall-type states have been predicted to occur in rotating Bose systems, and suggestions 
have been made how to manipulate and observe their fractional quasiparticle excitations. Due to a 
very rapid development in experimental techniques over the past years, experiments on BEC now 
appear to be close to reaching the quantum Hall regime. This paper reviews the theoretical and 
experimental work done to date in exploring quantum Hall physics in cold bosonic gases. Future 
perspectives are discussed briefly, in particular the idea of exploiting some of these strongly correlated 
states in the context of topological quantum computing. 
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I. INTRODUCTION 



Almost a decade ago, it was realized 



[\p\ that there exists an intimate theoretical connection between two seemingly 
very different physical systems: The quantum Hall effect(QHE) J], which occurs in a two-dimensional electron gas 
(2DEG) at very low temperatures and strong magnetic fields, and a rapidly rotating, dilute Bose-Einstcin condensate 3| 
of electrically neutral atoms. The key to this analogy is the observation that in two dimensions, rotation and a 
perpendicular magnetic field play a very similar role, making the two systems mathematically equivalent. This 
implies that, at sufficiently fast rotation, a Bose condensate is expected to enter a regime with strongly correlated 
states of the quantum Hall type, including quasiparticle excitations that obey fractional (anyon) statisticsjji. 

In recent years, the quantum Hall - rotating BEG analogy has been theoretically explored in great detail, starting 
with the prediction 5] of a bosonic Laughlin state[6|]. In addition, a large number of other incompressible states are, 



in principle, expected to occur, including the bosonic analog of the Jain sequencej?! and other Abelian states, but 
also non- Abelian states [s, 9|. 



with spin 



10| 



Moreover, a number of exotic states have been predicted in the case of rotating bosons 



Unfortunately, experiments have not yet reached this quantum Hall regime. However, there is reason to be opti- 
mistic, as the experimental development over the past decade has been astounding. The first experiments on rotating 
atomic Bose condensates were performed in the late nineties, and the first observation of a quantized vortex re- 
ported in Since then, Abrikosov lattices with hundreds of vortices have been produced, and present-day 



experiments 



13[ are close to the rotation speed at which this vortex lattice is predicted to melt and the system would 



enter the quantum Hall regime. The main obstacle is that at these rotation speeds, the system is close to the point 
where the centrifugal potential cancels the external harmonic trap and the atomic cloud would fly apart. There are, 
however, recent proposals how to get around this problem by modifying the confining potential, so that the quantum 
Hall regime may be reachable in the near future. 

In addition to being interesting in its own right, the prospect of producing quantum Hall- type states in cold atom 
systems may have long-term practical applications. One of the reasons for the recently revived interest in the anyonic 
excitations of the QHE is the theoretical proposal to use them in the context of quantum computing IJ] . This vision 



is certainly very far into the future. On the other hand, atomic systems may eventually turn out to be superior to 
the 2DEG quantum Hall system, as they allow for a very large degree of controlled tunability of various experimental 
parameters, including interaction strength and details of the confining potential. Moreover, these systems are well 
isolated and clean and thus less prone to decoherence than solid state realizations. 

This paper presents a (hopefully) comprehensive review of this, still active, fleld of research. Section ITTl gives a brief 
overview of the fractional quantum Hall effect. In Sec lIIII we give a general introduction to the subject of rotating 
Bose condensates, summarizing the experimental developments of the past decade and the theoretical understanding 



of how the system goes from a vortex lattice to the quantum Hall regime as rotation is increased. Sec lIVI explains 
the theoretical equivalence between a fast-rotating Bose gas and electrons in a strong magnetic field, along with 
some of the basic properties of the resulting many-body spectrum in the presence of interactions. An account of 
the literature on Abelian bosonic quantum Hall states is given in Sec|Vl most of this work is numerical or based on 
trial wave functions such as those of the composite fermion phenomenology [7[ and involves testing for the occurrence 
and stability of incompressible states (and their fractional excitations) at, e.g. the Jain fractions. This section also 
contains a brief discussion of the applicability of the composite fermion scheme for very low angular momentum 
states. Quite some work has been done to study the possible occurrence of non-Abelian quantum Hall states, which 
are particularly interesting in the context of topological quantum computing. This is accounted for in Sec lVIl which 
concludes that such states appear to be more prominent in rotating BEC than in the conventional QHE. Following 
this, we summarize various recent proposals how to design experiments capable of reaching the quantum Hall regime 
with present-day experimental techniques fSec lVIl] 1. Finally, we briefly discuss multicomponcnt Bose condensates 
fSec lVHl)) and round off with come concluding remarks and future perspectives in Sec lIXl 

II. THE FRACTIONAL QUANTUM HALL EFFECT - A BRIEF OVERVIEW 

The fractional quantum Hall effect (FQHE)^ is one of the most intriguing and most studied phenonema in con- 
densed matter physics during the past 2-3 decades. It occurs in two-dimensional, high-mobility electron systems 
(typically formed at the interface between two semiconductor crystals, e.g. in GaAs heterostructures) subjected to a 
strong magnetic field and low temperatures (in the millikelvin regime in present-day experiments). Figure [T] shows a 
sketch of a typical Hall experiment: A current Ix is passed through the sample along the x-direction and the resulting 
transverse voltage Vy measured for varying values of the magnetic field. Roughly speaking, the occurrence of a trans- 
verse voltage can be understood as being due to the deflection of the charge carriers in the presence of the external 
magnetic field, causing a build-up, or imbalance, of charge along the edges of the sample. In a purely classical picture, 
this leads to a linear relation for the Hall resistance R^y — Vyjlx as function of the magnetic field B. In the quantum 
Hall effect, however, the Hall resistance is quantized, 

^"^"^ Ix ~ ly e^' ^ ' 

where h is Planck's constant and e the electron charge. The number v takes integer values (integer quantum Hall 
effect) or equals rational fractions (fractional quantum Hall effect), and each allowed value of Rxy remains constant 
for a finite range of the magnetic field, as indicated in figure [51 At the same time, the longitudinal resistivity pxx 
equals zero, except for the transition regions between neighbouring plateaux. A third characteristic property of each 
of these quantum Hall states is that the system is incompressible, i.e. there is a gap between the ground state and 
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FIG. 1: Sketch 
passed through 
magnetic field. 



of the Hall experiment. The 2DEG is exposed to a strong perpendicular magnetic field B^. A current is 
the sample along the a;-direction, and the resulting transverse voltage Vy measured for varying values of the 
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FIG. 2: Sketch of the (integer) quantum Hall effect. The Hall resistance as function of the magnetic field is quantized, i.e. 
exhibits plateaux. (The result predicted by the classical Hall effect corresponds to a straight line through the centres of the 
plateaux.) The longitudinal resistivity is zero except at transitions between plateaux. Courtesy of D.R. Leadley, Warwick 
University 1997. 



15| : two years later, using even 
16| at v = 1/3. Since then, with 



17j. The quantization 



(bulk) excitations. The integer effect was discovered in 1980 by von Klitzing et al 
cleaner samples, Tsui and collaborators reported the discovery of the fractional effect 
the fabrication of ever-higher mobility samples, a large number of fractions have been observed 
of the Hall resistance turned out to be extremely exact (to at least ten parts in a billion), which has led to the 
introduction of a new standard of resistance, with the so-called von Klitzing constant Rk — h/e'^, roughly equal to 
25812.8 ohms, as the fundamental unit. 

Physically, the number in ([1]) corresponds to the Landau level filling fraction at the center of the corresponding 
plateau In other words, the IQHE occurs when an integer number of Landau levels is filled, while the FQHE 



is seen at fractional filling. As a consequence of this, the integer effect may be qualitatively understood in terms of 
a non-interacting electron picture 2]. The fractional effect, on the other hand, is a much more subtle phenonemon, 
taking place (mostly) in the lowest Landau level (LLL) and possible to understand only when the interactions among 
the electrons are taken into account. The theoretical explanation of the most prominent fractional states at i/ = 1/m, 
where m is an odd integer, was given by Laughlin in 1983, when he proposed his famous many-body wave function 
for the ground state ^ , 

eB 



Here, Zi = (xi +iyi)y j— are two-dimensional complex coordinates denoting the positions of the particles in the plane. 
It makes intuitive sense that this wave function does a good job in minimizing the Coulomb energy of the strongly 
correlated electrons in the plane - it contains the mth power of a Jastrow factor, which is a product of the 'distances' 
(relative coordinates) between all pairs of particles. The Jastrow factor approaches zero when any two particles try 
to come close to one another and thus, in a sense, helps to keep the particles apart jl^. This is a useful picture to 
keep in mind, as this factor will show up in various contexts later in the article. 

The Laughlin wave function describes a novel, inherently quantum mechanical state of matter, an incompressible 
quantum fluid. One of the most exotic properties of this state, and all other fractional quantum Hall states, is that 
it supports fractional excitations; Laughlin showed that its fundamental, charged quasiparticle excitations carry a 
fraction (1/m) of the electron charge and obey fractional (any on) statistics J|. The latter means that they are neither 
bosons nor fermions; when two such quasiparticles are exchanged in a counterclockwise manner, their wavefunction 
picks up a phase 

V'(r2,ri)=e""V'(ri,r2) (3) 

with a = l/m for the v ~ 1/m state. (Bosons and fermions would correspond to a = and 1, respectively.) Laughlin's 
theoretical explanation of the FQHE earned him the 1998 Nobel prize in physics, together with the experimental 
discoverers Tsui and Stormer. 

As mentioned previously, a large number of FQH states have been observed since the discovery of the = 1/3 
plateau. Most of these occur at odd denominator fractions, and many (but not all) belong to the "Jain sequences" 
u — n/{2np± 1) with n and p integers. There have been two main theoretical approaches to these states. In the 



Haldane-Halperin hierarchy picture 



20 



2l( . a QH state can give rise to a sequence of 'daughter states' as successive 
condensates of quasielectrons and/or quasiholes. The basic idea is that, once the system is in a Laughlin state and 
a sufficiently large number of quasielectrons or -holes have been generated (typically by changing the magnetic field 
away from its value at the center of the quantum Hall plateau), these quasiparticles themselves may form a strongly 
correlated state, in much the same way as the electrons form the Laughlin state. The result is a new incompressible QH 



ground state at a different filling fraction, whose quasiparticles may again condense to form the next 'daughter' etc. 

n 

The other approach is based on Jain's phenomenology of composite fermions 7]. The main idea of this construction 
is, roughly speaking, to attach an even number of vortices to each electron. These vortices effectively cancel a part 
of the external magnetic field, thus mapping the electrons into weakly interacting composite fermions which can then 
be thought of as moving in a reduced magnetic field. This picture provides a method to construct explicit trial 
many-body wave functions for the ground states at i' — n/{2np± 1), as well as their quasiparticle excitations. This 
approach has proven highly successful, producing wave functions with very high overlaps with the corresponding exact 



3. 



based on the use of conformal field theory 



states; we will get back to the details in section IV Al Recent work ; 
methods to construct hierarchical FQHE wave functions, illustrates that these two seemingly competing approaches 
are, in fact, very closely related Igs], I2J1. 

While almost all polarized FQH states observed to date occur at odd-denominator fractions, there is one known 
exception, namely, the gapped state at ly ^ 5/2. It is believed to be described by the so-called Pfaffian wave 
function proposed by Moore and Read^Sj. The Pfaffian is the exact ground state of a three-body repulsive interaction 
and describes a paired state very similar to a p-wave superconductor it is even more exotic than the states 
discussed above, in that its quasiparticle excitations obey non-Abelian fractional statistics 1^. This generalization of 
'conventional' (Abelian) anyons requires a degenerate set of d states with quasiparticles at fixed positions ri, r2, • • • 
r„ , such that an interchange of two quasiparticles i and j corresponds to a unitary operation in the subspace of these 
degenerate states, 

V-a - P^S^p. (4) 

Here, p*-*^^ is a d x d unitary matrix, and the set {ipa} denotes an orthonormal basis of the degenerate states. If 
the unitary matrices corresponding to different quasiparticle interchanges do not commute, the particles are said to 
obey non-Abelian statistics. The recent theoretical proposal that it might be possible to use such non-Abelian anyons 
in the context of topological quantum computing[14] has spurred great interest in the physics of the v = b/2 state. 
Moreover, there exist mathematical generalizations of the Pfaffian, the so-called Rezayi-Read (RR) or parafermion 
states [9|; although there are speculations that the recently observed QH plateau at v — 12/5 26] might correspond 
to a fc = 3 parafermion state, there is so far no unambiguous evidence of the existence of such, even more exotic, 
non-Abelian states in the quantum Hall system. 

Direct analogies of all the above (and many more) features of the FQHE are, in principle, expected to occur in rapidly 
rotating Bose gases and have been extensively studied in the literature in recent years. After a general introduction 
to rotating Bose condensates in the next section, the remainder of this article will be devoted to discussing these 
analogies in more detail. 
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FIG. 3: Image of an Abrikosov vortex lattice in a rotating BEG. From the JILA web page; courtesy of Eric Gornell. 



III. ROTATING BOSE CONDENSATES 

Bose-Einstein condensation of magnetically trapped alkali atoms was first achieved in 1995 '271, opening up many 
new directions of research on the border between atomic and condensed matter physics. Soon after these seminal 
experiments, people got interested in the rotational properties of atomic Bose condensates, and the occurrence of 



quantized vortices [28| was predicted 



The subsequent experimental development has been astonishing. The first 

11 1, and soon after by the Paris group [l^. 



ever vortex in such an atomic cloud was reported by the JILA group in 1999 
In the former, the vortex state was obtained by a direct imprinting of the 2t: phase shift onto the condensate, while 
the latter experiment used a mechanical stirring technique, with laser beams acting basically like a spoon in a cup 
of coffe elSCM . Following these experiments, the same stirring technique was used to create ever larger amounts of 
vortices 131. 3 



33|, which could be seen to organize themselves in triangular (Abrikosov) vortex lattices. A well- 
established technique to visualize, e.g. , such vortex arrays, is to perform absorption imaging along the rotation 
axis - a picture is taken after switching off the trap and allowing the cloud to expand for a fraction of a second. 
The vortices then appear as density dips in the image, as shown in figure [3] There is a third method which can be 
applied to further increase the angular momentum once the cloud is rotating; it is based on "evaporative spinup" , i.e. 
evaporating atoms with angular momentum smaller than average 
to create arrays with up to 200 vortices [; 



vortex cores 



36| 



M 



351 . Using this technique, it has been possible 
36j . and to study detailed properties such as the vortex modes Q, [s^] and 



An interesting question to ask, then, is what will eventually happen to the vortex array if one keeps increasing 
the angular momentum of the system and thus the density of vortices. One might expect that at some point, the 
vortex cores would start to overlap, as is the case in type II superconductors about to go normal. However, the 
picture that has emerged is quite different [ssjl . Since the confining potential in typical BEG experiments is harmonic, 
there is a limit where the centrifugal potential cancels the external potential and the cloud will get deconfincd and 
fly apart (see Sec lIVI for a more detailed discussion). When approaching this limit from below, the particles spread 



out, making the cloud more and more pancake-shaped, and the effective interparticle interaction becomes weaker due 
to the decrease in density. As we shall see in the next section, these are preconditions for the BEC to be equivalent 
to a system of particles in the lowest Landau level. The entrance into the LLL regime is signalled by a shrinking of 
the vortex cores starting around the rotation frequency at which the size of a vortex core becomes comparable to the 
spacing between vortices - the ratio between the size of a vortex core and the area occupied per vortex, saturates to a 
constant. This behaviour was predicted by Baym and Pethick , 



39( 1 and confirmed in recent experiments 13] where 



the transition to the LLL regime was observed to occur at around 98% of the deconfinement limit. Eventually, at 
even higher rotation, the vortex lattice is expected to melt. This melting transition has been studied theoretically by 



several groups 



4C 
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42| and is predicted to occur around ^ 6 — 10, where the filling factor v is the ratio of boson 
density to vortex density. Beyond this point, the system enters a regime of homogeneous, strongly correlated states 
of the same nature as those in the fractional quantum Hall effect. 

So far, experiments have not actually reached this quantum Hall regime. Rotation frequencies of more than 99% 
of the deconfinement limit have been achieved[l3], which is believed to be close to the vortex melting transition. The 
main practical problem is to push the rotation further without passing the point where the cloud flies apart. There are 
several very recent proposals of ways to avoid this problem, basically by modifying the external confinement. With 
these modifications, it may well be possible to reach the quantum Hall regime with presently available experimental 
techniques. We shall discuss these novel ideas in some more detail in Sec lVIII Meanwhile, the next three sections 
summarize the theoretical analogies between rapidly rotating Bose condensates and the quantum Hall effect, simply 
assuming that the system is in the lowest Landau level regime. 



IV. ROTATING BOSONS AS A LOWEST LANDAU LEVEL PROBLEM 

The basic insight, providing the analogy between rotating Bose condensates and the quantum Hall system, is that 
the Hamiltonian of a rotating system of harmonically confined, neutral particles is essentially equivalent to that of 
charged particles in an external magnetic field. We start by making this statement more precise and discussing under 
what circumstances a rotating Bose gas can be mapped to a lowest Landau level problem. The second part of this 
section presents some of the general properties of the corresponding many-body energy spectrum. 



A. Mapping to the LLL 



Let us consider a system of N spinless bosons with mass to in a harmonic trap of strength lu, rotating with angular 
frequency and interacting via a short-range (delta function) potential Hj. In a rotating frame the Hamiltonian can 



10 

be written as 



^ r -2 1 

2m 2 



ni^ + Hi (5) 



where Lz denotes the angular momentum around the rotation axis. (For simphcity, we wiU set h — 1 whenever there 
is no risk of confusion.) Separating out the planar (x, y) part and completing the square inside the brackets, one can 
rewrite Eq.® as 



' +{uj- n)L, + Hi (6) 



^ r 1 2 

■^-^ 2m V / 



with A = mLu{—y, x), \\ denoting the planar {x, y) part of the Hamiltonian, while Hho{z) denotes the perpendicular (z) 
part of the harmonic oscillator potential. This is how the formal link to the quantum Hall system comes about: We see 
that the planar part of H takes the form of particles moving in an effective "magnetic" field -Be// = V x ^4 = 2niujz. 
The quantum mechanical one-body spectrum of this part of the Hamiltonian is given by the so-called Landau levels 
(see, e.g. , ^), with energy = ('^ + 5) where n = 0, 1, 2, and ujc — 2ti;. Each Landau level is degenerate 
in angular momentum, the number of states per Landau level being equal to the number of (effective) flux quanta 
piercing the plane. The single-particle wave functions in the symmetric gauge chosen here, can be expressed as 

r/„,™ = ^„,™e-l^l'/4z™i™(y), (7) 

where n is the Landau level index, m denotes the angular momentum, Nn.m is a normalization factor, L™ are the 



associated Laguerre polynomials, and z = \/2muj{x + iy) is again a (dimensionless) complex coordinate denoting the 
particle position in the plane (note the change in notation as compared to ([5])). Now, the interaction is assumed to 
be weak in the sense that it does not mix different harmonic oscillator levels. We will be interested, for a given total 
angular momentum, only in the lowest lying many-body states (the "yrast" band). In this limit, the model may be 
rewritten as a lowest Landau level (LLL) problem in the effective "magnetic" field i?e// — "imuj (and — for the 
harmonic oscillator in the z-direction). The Hamiltonian then reduces to the form 

H = {LU-n)L + gJ2S^i^i-i^j) (8) 

i<j 

where L denotes the total angular momentum, L = li = L^. The single particle states spanning our Hilbert space 
(the lowest Landau level) are thus (omitting normalization factors) 

Vo,,n = z"^e-''/* (9) 

A general bosonic many-body wave function ^{zi, ...zn) with good angular momentum can thus be expressed as a 
homogeneous, symmetric polynomial in the coordinates {zi}, times the exponential factor exp{—'^^ |zip/4) (which 
will be suppressed throughout most of this paper for simplicity) ; the degree of the polynomial gives the total angular 
momentum of the state. 
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FIG. 4: Sketch of the single-particle Landau level spectrum. The slope of each level equals uj — fl. 

In the theoretical approach employed in the following sections, the system is simply assumed to be sufficiently dilute 
{i.e. weakly interacting) to be in the lowest Landau level at all angular momenta; the general strategy will be to look 
for the lowest (interaction) energy state within this subspace for given L and any fixed il. The experimental situation 
is somewhat different - we saw in Sec lIIII that in order to become sufficiently spread out to be in this dilute regime, 
the cloud, in present-day experiments, has to rotate faster than Q w 0.98 a;. Moreover, a more natural picture in an 
experimental setting is to think of as fixed, while the system selects the ground state angular momentum such as 
to minimize Ei + (uj — fl)L. Then, in order to avoid Landau level mixing, the number of particles and/or (w — fi) 
have to be small (see figH]). 



B. Yrast spectra 

A convenient way of studying the many-body properties of a rotating boson system is to display its yrast 
spectrum 4^ , where the lowest many-body energy eigenvalues are plotted as a function of total angular momen- 
tum. An example is shown in Fig. [SJ which was obtained by an exact diagonalization of the Hamiltonian ([5]) of the 
previous section for four particles, with the lowest Landau level restriction imposed. In the absence of interactions, 
the system is highly degenerate - the degeneracy for a given total angular momentum L corresponding to the num- 
ber of ways L quanta of angular momentum may be distributed among N bosons. This degeneracy is lifted by the 
short-range repulsion, leading to the energy band seen in the figure (the spectrum is shown for Q ^ u, i.e. shows 
purely the interaction energy, cfr. Eq.®). The line connecting the lowest states at different angular momenta is 
commonly called the yrast line. A number of basic properties of the system may be read out of this spectrum. First 
of all, one notices that the lowest possible energy decreases with increasing angular momentum, starting from L = 
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FIG. 5: Yrast spectrum, i.e. many-body interaction energy as function of total angular momentum, for four bosons in a 
harmonic trap with delta function repulsion. This figure was first published in 



where all bosons sit in the lowest angular momentum state; this is due to the particles' ability to spread out more 
in the plane, as more and more angular momentum states become available to them (or, in other words, because the 
particles interact only for zero relative angular momentum). In particular, one notices that at L = N(N — 1) = 12 
and above, the ground state has zero interaction energy. This is the first immediate consequence of the analogy to 
quantum Hall physics: Just as in the quantum Hall effect, it is possible to construct particularly well correlated states 
of the Jastrow form, i.e. 



(10) 



i<j 



where m is a positive integer, and /({zi}) is some homogeneous, symmetric polynomial in the N coordinates {zi}. 
Since the repulsive interaction between the bosons is zero-range, any state of this form will have zero interaction 
energy. Since the total power of a Jastrow factor Y[i<j i^i ~ ^j) equals the number of pairs of particles, N{N — l)/2, 
the Jastrow part of the wave function pop contributes an angular momentum Lq = mN{N — 1). Therefore, the 
smallest angular momentum at which a state of the type (jlOp can exist, is N(N — 1). At this angular momentum, 
the exact, non-degenerate ground state of the system is given by the Bose-Laughlin state 



(11) 



n 



i<3 



as was first pointed out in Ref. [5]. The degeneracy of the zero-energy yrast states dX L > Lq corresponds to the 
number of ways L ~ Lq angular momentum quanta can be distributed among the N particles; alternatively, this 
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degeneracy can be found 
level 



45| by exploiting the fact that the wave functions (fTO|) describe anyons in the lowest Landau 
46| | obeying Haldane's exclusion statistics |47j with statistics parameter g ~ 2. 
Another important feature of the yrast spectrum is that for each state, there is a set of 'daughter states' at all higher 
L, with the same (interaction) energy. These are simply center-of-mass excitations of their 'parent' state Lsl- At a 
given L, all states corresponding to center-of-mass excitjitions of lower-L eigenstates, are orthogonal to the subspace 
of "new" states. According to Trugman and Kivelson[48l|. the latter subspace consists of translation invariant (TI) 
polynomials, i.e. polynomials which arc invariant under a simultaneous, constant shift of all coordinates, Zi — > + a. 
A convenient basis for the TI subspace are the elementary symmetric polynomials of degree L, SL{{zi]) (2 < i < N), 
defined as 



sl(zi, . . . , zat) = 5 {ziZ2...zl} . 



(12) 



with z = z — Z where Z = ZijN is the center of mass, and S denotes synimetrization of the product over all N 
particle coordinates. Together with si({zi}) = zi -f Z2 + ••■ + -^w = NZ., the elementary polynomials (|12p thus span 
the entire space of symmetric, homogeneous polynomials in the lowest Landau level. We will make explicit use of the 
basis (I12|) when discussing low-angular momentum states in Sec lVDl 

The yrast spectra of rotating Bose condensates in the weak interaction (lowest Landau level) regime have been 
studied extensivel y in the literature, using a variety of methods including analytical studies 
Pitaevskii) theory 



49[, and exact numerical diagonalization 
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44|, mean field (Gross- 



. For example, Reimann et al[5l!| demonstrated 



how the presence of localized vortices in a rotating boson cloud is revealed by periodic cusps in the yrast line of the 
exact many-body spectrum. However, the dominating line of attack in studying the analogies to quantum Hall physics 
has been the use of various types of trial wave functions. These include bosonic versions of Laughlin- and Jain- type 



wave functions 
wave functions 
sections. 



and variations thereof 
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53[, as well as the bosonic counterparts of non-Abelian quantum Hall 



8|, |9|. The results of these trial wave function studies are reviewed in more detail in the following 



V. ABELIAN QUANTUM HALL STATES 

Most known fractional quantum Hall states are Abelian. These include the Laughlin- and Jain states and more 
generally all states which, in the hierarchy picture, are generated by any sequence of quasielectron- and quasihole 
condensates, as was discussed in section |lTl The bosonic counterparts of the Laughlin- and Jain states have been 
extensively studied in the literature, mainly numerically. Apart from exact diagonalization, a particularly widely used 
technique is the composite fermion approach; this line of work is reviewed in the first and main part of this section. 
In addition, we briefly discuss the anyonic quasiparticlcs of these states, the bosonic hierarchy, and the possibility to 
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apply the composite fermion formalism at very low angular momenta. 



A. The Jain sequence and composite fermions 

We argued in the previous section that translation invariant (TI) states play a special role in the rotation spectra. 
A particularly useful scheme of constructing TI trial wave functions that has been widely exploited in quantum Hall 
physics, comes from the phenomenology of composite fermions (CF) 5^. Composite fermions were first introduced 
by Jain[7i] and have proven very successful in describing FQH states, quantum dots in high magnetic fields [55] and, 
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56 



57| . In quantum Hall physics, the 



as we shall now discuss, highly rotational states of Bose condensates [5, 
basic picture of Jain's construction is, roughly speaking, that an even number of vortices is bound to each electron. 
Each of these vortices effectively cancels one flux quantum of the external magnetic field, and the electrons are thus 
mapped into weakly interacting composite fermions, which can be thought of as moving in a reduced magnetic field. 
Technically, "attaching a vortex" means multiplying the wave function by a Jastrow factor, 

l[{z, - z,). (13) 

We already mentioned that the Jastrow factor has the effect of keeping the particles apart - it goes to zero if any 
two coordinates Zi and Zj approach each other. Therefore, it takes care of much of the repulsive interaction between 
the particles. In the simplest approximation, the so-called non-interacting composite fermion (NICE) approach, the 
composite fermions are thus simply assumed to be non-interacting. This kind of considerations led Jain to construct 
trial wave functions a Slater determinant of (non-interacting) composite fermions in the reduced magnetic field, times 
an even power of Jastrow factors. In the case of bosons, whose wave function has to be symmetric rather than 
antisymmetric, the construction is modified by instead binding an odd number of vortices, mapping the bosons to 
weakly interacting composite fermions. In other words, bosonic trial wave functions with angular momentum L are 
constructed as non-interacting fermionic wave functions with angular momcntmn L —pN{N— l)/2, multiplied by an 
odd number p of Jastrow factors, and projected onto the LLL, 

^L=V ^fs{z^, z,) n(z, - j . (14) 

Here, fs denotes a Slater determinant consisting of singlc-particlc wave functions ([7]). The LLL projection V amounts 
to the replacement Zi —^ 2d/dzi in the polynomial part of the wave function ~ the recipe is to replace all z:s with 
derivatives in the final polynomial, after multiplying out the Slater determinant and the Jastrow factors and moving 
all z:s to the left. It has been shown 7] that with this projection method, the single-particle wave functions in the CE 
Slater determinant may be written as 

r?„i = I > -n (15) 
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with all derivatives acting only to the right. As this method tends to get computationally heavy in numerical 
calculations with many particles and a large number of derivatives, somewhat different methods of obtaining LLL 
wavefunctions have been employed in most of the CF literature 7]. These, too, are often referred to as projection. 
Nevertheless, in this paper, "projection" will refer to the above "brute force" procedure. 

Before summarizing the results obtained in the literature, let us illustrate the method on two simple and well-known 
examples in the QH regime: First, consider the case L — N{N — 1). Taking p = 1, the Slater determinant fs has to 
contribute an angular momentum N{N — l)/2 and is given by putting all CFs into the lowest CF Landau level, from 
I = Q to I = N - I, 
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We immediately see from Eq. (jl4l) that the full wave function is simply the bosonic Laughlin wave function (jll|) 
with angular momentum L = N{N — 1). Next, consider the angular momentum N{N — 1) ~ N, corresponding to 
a "quasiatom" (the bosonic counterpart of a quasielectron) at the center. While a quasihole can be seen as a local 
depletion of the quantum Hall liquid, a "quasiatom" corresponds to a local contraction, with fractional surplus charge 
(or rather particle number in our case of neutral bosons) 1/2. In the CF language, a trial wave function for such a 
quasiparticle excitation at the center, i.e. with minimum angular momentum, is obtained by exciting one composite 
fermion to the second CF Landau level, leading to the Slater determinant 
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We see that it is the excited composite fermion, with single particle wave function ~ z (to be replaced by a derivative 
upon LLL projection) that causes the reduction of angular momentum as compared to the Laughlin state. One obtains 
the full trial wave function (again, apart from the exponential factor) 
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with i/jl denoting the Laughlin state (fTTj) . This wave function has very high overlap with the exact one (e.g. , 99.7 



% for 4 bosons 



45|). Its fermionic counterpart has been proven to capture correctly both the fractional charge and 
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59|, and the same can be expected to be the case for this bosonic 



the anyonic statistics of the QH quasielectron 
version. 

Trial wave functions for other yrast states are constructed in similar ways. The lower the angular momentum, 
the larger the number of derivatives. Had we filled up the second CF Landau level, with equally many composite 
fermions as the first, we would have obtained a trial wave function for a new incompressible ground state, with filling 
fraction 2/3. In general, ground states of the principal Jain sequence v = n/{n + 1) - the bosonic counterpart of 
the well-known principal Jain sequence v — n/ {2n + 1) in the FQHE - are described as v* ~ n integer quantum 
Hall states oi p — 1 composite fermions, while quasihole- and "quasiatom" excitations are described by removing a 
CF from a filled CF Landau level and adding a CF to an otherwise empty CF Landau level, respectively. Of course 
this construction can be generalized in the usual way by attaching a larger (odd) number p = 3, 5, ... of vortices to 
each boson. Note however, that for a pure delta function interaction, the ground states at the corresponding filling 
fractions v = n/ {np+ 1) belong to the highly degenerate subspace of zero-energy states discussed in Sec II V Bl making 
this construction less relevant. Adding higher derivatives of the delta function to the Hamiltonian (or, equivalently, 
using Haldane's pseudopotentials[2]) lifts this degeneracy, and the CF construction with p > 1 again provides good 
wave functions for the resulting ground states [4^. 

The ideajjf applying the CF phenomenology in the context of rotating Bose gases was first tested by Cooper 
and Wilkin[^] and by Viefers et al[|45!] in disk geometry 2 1 for small numbers of particles, by comparing to exact 
diagonalization results. It was shown that the approach reproduces many prominent features of the yrast spectrum, 
such as the locations of cusps in the yrast line; overlaps with the exact solutions for a number of yrast states were 
computed for up to ten par ticles and shown to be large - typically 99% for five particles. Later, several more systematic 
studies were performed 



56 



57| in spherical geometry 



The advantage of this theoretical approach, in which the 
particles move on the surface of a sphere with a radial magnetic field produced by a magnetic monopole at the center, 
is that the sphere has no boundaries. While edge effects play an important role for small systems in the plane, this 
geometry thus allows for the 'simulation' of homogeneous bulk states even for the relatively modest particle numbers 
accessible to numerical calculations. Instead of the complex coordinates Zi discussed so far, the particle positions on 
the sphere are parametrized by the polar angles {9i, 4>i), or more conveniently, by the spinor coordinates 

u, = cos(6',/2)e*"^'/2^ = sm{^^/2)e-"''^/^ . (19) 

For example, the Bose-Laughlin wave function pT|) takes the form 

^ IT ~ ^j'^if ' (20) 
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and other wave functions may be translated from the plane to the sphere in a similar way. Performing exact di- 
agonalizations for up to 12 particles on the sphere, Regnault and Jolicoeur 5^ found evidence of the occurrence 
of incompressible (gapped) states at the principal Jain fractions v = ^, f , f , f , as well as excited states in general 
agreement with the CF phenomenology. Moreover, going away from pure delta function interaction by adding a 
higher-order pseudopotential (as discussed above), they found evidence of an incompressible state aX v — 2/5. This 
state is not part of the principal Jain sequence; rather, it is the bosonic counterpart of the 2/7-state in the FQHE. 
Provided it is Abelian, it can be understood, in the hierarchy picture, as resulting from a condensate of quasiholcs in 
the V = 1/2 Bose-Laughlin state; in the CF picture, it belongs to the negative p — ?> Jain sequence, v — n/{3n — 1), 
with n = 2. Although the interaction in typical experiments is dominated by s-wave scattering, there exist methods 
to introduce and enhance an additional d-wave interaction 



6C 



61 



62l |. This may provide a possibility, at least in 



principle, to observe the 2/5 state (and other states with i/ < 1/2). Alternatively, one might use a system of atoms 



with permanent dipolar interaction 
Chang et al 



63[ such as chromium 64]. Following up on the work of Regnault and Jolicoeur, 



57[ performed a direct comparison between exact diagonalization results and those predicted from the 



CF approach, computing energies as well as overlaps between exact and CF wave functions for the ground states and 
low-lying excitations at v = i, | and |. They found that the non-interacting composite fcrmion approach correctly 
predicts the incompressibility of the ground states at = i , | and | and produces excellent overlaps (over 97% for 
up to 10 particles) for the ground state and excitations at = i as well as the ground state at = |. However, for 
increasing n, the NICF approximation gets progressively worse, for short-range as well as Coulomb interaction, pro- 
ducing considerably poorer overlaps than in the principal Jain sequence of the electronic FQHE. In the CF language, 
the interpretation is that 'residual interactions' between the composite fermions play an important role. In particular, 
in the limit n — > oo, i.e. v ^ \ (the bosonic counterpart of the metallic v — ^ state in the quantum Hall system), 
the ground state of the system can not be described by a Fermi sea but rather appears to be a non- Abelian state, the 
bosonic version of the Moore-Read PfafSan ^. We shall get back to this point in Sec lVIl 

Additional evidence for the strongly correlated nature of the states at = 1/2, 2/3 and 3/4 was given by Cazalilla 
et al[65[ who studied the low-energy edge excitations of harmonically confined, rapidly rotating few-boson systems. 



According to Wen 6^, the 'topological order' of a bulk quantum Hall state, implying its filling fraction, as well as the 
charge and statistics of its quasiparticle excitations, is reflected in the properties of its edge excitations. Performing 
exact diagonalization studies for up to seven particles, Cazalilla et al showed that the number of edge modes is 
consistent with that predicted by Wen's theory, for the states at v — 1/2, 2/3 and 3/4. 
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B. Anyonic excitations 



A particularly interesting aspect of the fractional quantum Hall effect is the existence of fractionally charged 
quasielectron- and quasihole excitations ^ , which are expected to obey fractional (anyonic) statistics 4 1 . Obviously, 
the same type of quasiparticles should occur in the bosonic quantum Hall system, with 'charge' replaced by particle 
number. In the simplest case of the v = 1/2 Laughlin state, the quasihole would be a vortex with local lack of density 
corresponding to half an atom, and would obey semionic statistics, i.e. quantum statistics 'halfway' between bosons 



and fermions. Such a quasihole, located at zq, is described by the wave function 

i^qh{{Zi\) = W{Zi - ^o) W{Zk - Zlf . 

Paredes et al 



(21) 



k<l 



671 suggested that, in principle, such quasiholes can be created by piercing the Bose-Laughlin state 
locally with lasers. Adiabatically moving such a laser would then "drag" the quasihole along, enabling controlled 
interchange of pairs of quasiholes and thus a direct measurement of the anyonic phase tt/2 picked up under exchange. 



The latter would be particularly interesting - despite very promising recent experimental progress 
FQH system, a direct and unambiguous measurement of fractional statistics is still lacking. 



in the electronic 



C. Hierarchy 



In addition to the Jain states discussed above, the analogy with the FQHE in 
of (Abelian) hierarchical states that do not belong to the principal Jain sequence 



3, 



pri ncip le 
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predicts a large number 



69|. An example of such a 



state in the 2DEG is the one recently observed at = 4/ll[70|, whose bosonic counterpart would he v — 4/7. Trial 
wave functions for these hicrachy states, or at least those corresponding to "quasiatom" (as opposed to quasihole) 
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68 



7l|; although this construction is 



condensates, may be constructed using conformal field theory techniques 
well-defined, its validity will eventually have to be determined by numerical tests of the resulting wave functions. As 
discussed above, even the principal Jain states do not describe the bosonic system as accurately as is the case in the 
2DEG, and the same may be the case for the general hierarchy construction. Moreover, as we have seen, trial wave 
functions for ground states at < 1/2 are of interest only in systems with scattering in higher partial waves. 



D. Low angular momenta — a digression 



Let us briefly address some interesting analytic results in a case that is very far away from the quantum Hall 
regime, namely the very lowest angular momenta up to the single vortex, 2 < L < N . Within the lowest Landau 
level approximation, exact ground state wave functions for all angular momentum states in this interval were derived 



some years ago 



73|. They are given by the elementary symmetric polynomials s^^Zi) where Zi 



Z and 
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FIG. 6; Overlap between the CF trial wave function sd, L — N and the exact one given in Ea. (|22|l . as a function of the number 
of particles. Note the range of numbers on the y-axis - the overlap equals 99.5% already for 10 particles. 

Z — J2i Zi/N is the center-of-mass coordinate, 

E {zp^~Z){zp,-Z)---izp,-Z). (22) 

pi <P2<---<PL 

For example, 'ipL=2 — S [{zi — Z){z2 — Z)], with S denoting symmetrization over aU particle coordinates; for L = N 
this expression reduces to ipL=N = Yiii^i ^ ■^)- Since in present-day experiments, the lowest Landau level ap- 
proximation is certainly not valid at these lowest rotational states, the results discussed in this subsection may be 
somewhat academic. However, there is an interesting connection between the exact wave functions (j22p and those 
following from a naive application of the composite fermion construction. A priori, one would certainly expect the 
CF construction to fail in this regime, at least if the usual qualitative picture of composite fermions were to be taken 
literally. According to this picture, it is the flux attachment (i.e. the factor Yii^i ~ ^jY) which makes the composite 
fermions weakly interacting, justifying the NICF approach. Note, however, that since the Jastrow factor itself has an 
angular momentum of mN{N — 1), the construction of CF trial wave functions at L ^ N involves applying 0{N'^) 
derivatives 
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74l | . One would expect that these derivatives acting on the Jastrow factor destroy most of the good 
correlations which are at the very heart of the CF construction. It is therefore intriguing that in this regime, a naive 
application of the NICF scheme produces wave functions whose overlap with the exact ground states (|^^ is not only 



74| the single vortex state, L = N, was studied in detail. 



large, but increases with increasing particle number. In Ref. 
and numerical calculations for up to 43 particles showed that the overlap between the CF trial wave function and the 
exact analytical result (1^^ approaches unity for large N, with the difference decreasing as ~ 1/A^ (see FiglB]). Further 
analysis 75j showed that this is not an artifact of the L — N state. Numerical tests up to iV = 7 for L — N — 1 and 



20 

up to A'^ = 8 for L = N —2 show the same tendency, i.e. . overlaps increasing with particle number. (For fixed N , on 
the other hand, overlaps tend to decrease as L decreases.) In fact, the CF wave functions can be shown to have an 
analytical structure strikingly similar to the exact ones. The simplest example is the CF state at i = iV = 4, which 
can be expressed as 

^P^^M=, - E n - Z^'') , (23) 
k=l 1=1 

with Z'^^'^ denoting the 'incomplete' center-of-mass coordinate '^j^kZj/{N — 1). The situation becomes more com- 
plicated for higher N and lower L, with more and more coordinates 'missing' in the center of mass (making the task 
of analytically proving the numerical results highly non- trivial), but apart from this, the general structure (|22p is 
reproduced by the CF construction. 



VI. NON-ABELIAN QUANTUM HALL STATES 

Although the zoo of Abelian states discussed in the previous section is very rich, it does not exhaust all possibilities 
of gapped quantum Hall states. The possible existence of non- Abelian states in the quantum Hall system was already 
pointed out in section |TT1 in this section we shall discuss the bosonic analogies of these states. Non- Abelian quantum 
Hall states have received great attention lately, due to the recent proposal to use their quasiparticle excitations in the 
context of topological quantum computation 14]. The main advantage of this scheme is its intrinsic fault tolerance 
- quantum information stored in states with multiple non- Abelian quasiparticles is 'topologically protected', i.e. 
immune to local perturbations. To actually perform a quantum computation would involve creating a state with 
a given number of quasiparticles at certain positions and performing controlled braiding operations by physically 
dragging quasiparticles around one another in a specified manner. One may thus speculate that, eventually, rotating 
Bose condensates may provide better candidates for topological quantum computing than the conventional QHE: 
In general, cold atomic systems are much easier to control and manipulate, with a high tunability of experimental 
parameters, and as already mentioned, there exist theoretical proposals 67| of ways to drag quasiparticles through 



the condensate. Moreover, non- Abelian states may actually be more prominent in the bosonic case; numerical studies 
indicate that both the Pfaffian and other parafermion states can be expected to occur already in the lowest Landau 
level. So it is obviously of interest to study the possible occurrence of such states in the bosonic system. The next 
two subsections summarize the work that has been done in this direction. 



A. The Pfaffian &t v = 1 



As was pointed out in Sec IV Al v ~ \ is, the bosonic counterpart of the half-filled Landau level in the electronic 
FQHE. In the latter, there is no quantum Hall plateau around v = 1/2 [2! [t^; rather, the system displays a compressible 
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state, behaving like a degenerate gas of (almost) free electrons in zero external magnetic field. This behaviour has 
been explained [t^I in terms of the non- interacting composite fermion model - at j/ = 1/2 there are exactly two 
external flux quanta per electron, so after binding two flux quanta each, the composite fermions are left in zero 
effective field and form a Fermi sea. On the other hand, at v = 5/2, one does find an incompressible state which is 
believed to be described by the Moore-Read wave function, the quantum Hall analog of the paired state in a p-wave 
superconductor Q|; the picture is that the presence of the two filled lowest spin-subbands effectively modifies the 
interaction between the electrons in the half-filled topmost Landau level, leading to pairing. 

In the bosonic system, the situation is qualitatively different, in the sense that there appears to be no compressible 
state at J/ = 1. Rather, there is quite substantial numerical evidence that the ground state corresponds to the bosonic 
version of the Moore- Read state, 

where the Pfaffian is defined as 



Zi Zj 



(24) 



Pf 



1 



A 



1 



1 



1 



(25) 



(zi - Z2) (z3 - Z4) {zn-i-zm)_ 
with A denoting antisymmetrization over all coordinates. The possibility of a Pfaffian state at = 1 was first suggested 



by Wilkin and Gunn 



78[; work by the same group 



40| later showed the existence of an incompressible groundstate 
dX V = 1, and reported large overlaps (more than 96% for six particles) between this state and the Pfaffian trial 
wave function; overlap calculations by Chang et a llSTII for up to 16 particles confirm this picture. Further numerical 
evidence was given by Regnault and Jolicoeur 56|, |79| , both for the ground state (including evidence of pairing from 
the form of the two-particle correlation function) and for quasihole excitations (correct degeneracies and high wave 
function overlaps). 

It appears clear that the existence of the Pfaffian state at v = 1 \s not an artifact of the short-range interaction - 
on the contrary, introducing a Coulomb interaction between the bosons even increases the overlap[57]. (Though the 
introduction of a strong d-wave component in the interaction may destroy the state[56]). In other words, the fact 
that the bosons have a stronger tendency of pairing than their fermionic counterparts, appears to be mainly due to 



their quantum statistics. 



B. Parafermion states 



An important difference between the system at hand and the 2DEG is that the bosonic system allows for states 
with v > 1 that are entirely in the lowest Landau level, due to the absense of Pauli blocking. It is thus of interest 
to understand what happens in the interval up to ~ 6 where the system is expected to enter the Abrikosov vortex 
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lattice regime [40j|. It was first suggested by Cooper et al[40j that in this interval, at v = k/2; k = 3,4,5, ... one may 
find a sequence of non-Abelian incompressible states described by the parafermion wave functions introduced by Read 

n 



and Rezayi|2]- They can be represented as 



N/k N/k 

n (^^-^.)' n (^'^■~^')' 

i<jeA k<l£B 



(26) 



where the system is divided into k groups {A,B, ...) each containing N/k particles, and S denotes symmetrization 
over all coordinates. The Laughlin state is recovered as the special case fc = 1, while the expression for fc = 2 is an 
equivalent way of writing the Pfaffian (|24p . Generalizing the pairing in the Moore-Read state, these wave functions 
describe states with fc-particle clustering; they are the exact zero-energy eigenstates of a {k + l)-particle delta function 
interaction. Performing numerical calculations on the torus. Cooper et al found large overlaps between the Read- 
Rezayi states ((26)) and the ground states at = fc/2 for 1 < k < 6 and also recovered the correct ground state 



degeneracies. Regnault and Jolicoeur 



56 



79| later took these calculations to larger systems (in spherical geometry) 
to see if this picture continues to hold as one approaches the thermodynamic limit. Their results remained somewhat 
inconclusive but indicated that for a pure delta function interaction, overlaps quickly decrease as k and the number 



8l| that the = 3/2 parafermion state is 



of particles are increased. On the other hand it was shown by Rezayi et al 
stabilized by introducing a moderate amount of longer-range interaction; similar conclusions were reached for v ^ 2 in 
a very recent paper by Cooper and Rezayil82']. In principle this may be achieved in a system with dipolar interactions 
or a moderate rf-wave component. What makes these states particularly interesting is that A; = 3 is the smallest /c- value 
among the parafermion states for which the non-Abelian statistics support universal quantum computation [l^ 



VII. BEYOND HARMONIC POTENTIALS - WAYS TO AVOID THE DECONFINEMENT PROBLEM? 



We have seen that one of the main practical obstacles to actually reaching the quantum Hall regime experimentally, 
is that, for the usual harmonic confinement, rotation speeds exceedingly close to the deconfinement limit are required. 



3, 



The present record, with rotation a,t Q > 0.99a;[13j, while having reached the lowest Landau level, still lies clearly 
within the Abrikosov lattice regime. An obvious way to be able to rotate the cloud faster than uj without it flying 
apart, is to modify the confining potential. This section summarizes a few such proposals. 



A. Quartic potentials 



Several theoretical studies have addressed the effect of adding a small quartic term to the trap 38[ . Among the 
predictions for the vortex array regime, are the occurrence of singly quantized vortex arrays with a hole in the middle 
or, at very high rotation, a single, multiply quantized vortex at the center of the trap. In experiment, such an 
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anharmonic trap has been created[30|, l8j| by applying a blue detuned laser propagating along the axis of the trap. 
This effectively amounts to adding a Gaussian potential U{r) ^ ?7o exp(— ar^) where r is the (planar) distance from 
the axis of rotation, and the constants Uq and a are given by the parameters of the laser. For small ar'^ this potential 
is well approximated by quadratic (giving a small correction to the original harmonic trap) + quartic. 

In a very recent numerical study, Morris and Feder 8J] propose that using this type of quartic potential would make 
it possible to attain the Bose-Laughlin state (and other quantum Hall states) with presently accessible rotation rates. 
They show that the inclusion of such a potential tends to lower the critical rotation frequencies at which the quantum 
Hall states are expected to occur. Moreover, they predict that fine-tuning of the Gaussian parameters (depending 
on particle number and interaction strength) is necessary in order to avoid destroying the Laughlin state, but that 



85j. In particular, the required experimental 



the required values of these parameters are within experimental reach [; 
parameters should become more easily accessible if the number of particles in the condensate is reduced; this may 
be achieved 86[ by adding a ID optical lattice along the axis of rotation, splitting the condensate into an array of 



independent quasi-2D BECs. 



B. Optical lattices 



831 ■ They 



A somewhat different, very recent theoretical proposal involving optical lattices, is due to Bhat et al 
suggest to include a co-rotating optical lattice (in the tight binding regime) in addition to the harmonic potential, 
keeping the system confined even at critical rotation velocity Q — to. In addition to avoiding deconfinement of the 
atom cloud, this model displays intriguing physical properties: Mapping the system to a Bose-Hubbard model, the 
authors show that the rotation introduces phase factors in the effective hopping term, Hhop ^ o|aje~"^'^ -I- h.c, 
where the phase depends on the rotation velocity and particle mass as well as the lattice spacing. The linear response 
of the system to a potential gradient (tilt of the lattice) shows quantum Hall-like features even for a single particle 
(and similarly for two particles). The authors give the following, qualitative explanation of the analogy to quantum 
Hall physics: The lattice, with its tunneling barriers, in some mean field sense mimics the repulsion experienced by a 
single particle from the rest of the 2DEG, with the inaccessible regions (maxima) of the lattice corresponding to the 
positions of the other electrons. Moreover, the phases picked up by a single particle when moving around the lattice, 
simulate the effect of the correlation holes (or vortices) at these 'electron positions'. But clearly a lattice potential 
cannot support a liquid state, so the exact correspondence between this system and the FQHE remains to be fully 
clarified by further studies. 

Finally, it is worth mentioning that there have been other theoretical proposals to create a quantum Hall effect 
for bosonic atoms, involving non-rotating optical lattices [88[, were the magnetic field is simulated e.g. by means of 



laser-induced hopping[89| . A particular advantage of optical lattices is the extremely large degree of controllability. 
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not only of the amount of flux per lattice cell but also the amount of disorder in the system. 

All the proposals discussed here, claim that the relevant model parameters are more or less within reach of present- 
day experimental techniques. Given these latest ideas, one can certainly hope for exciting experimental developments 
in the near future! 



VIII. MULTICOMPONENT BOSE CONDENSATES 



So far we have restricted ourselves to single-component condensates of spinless (or scalar) bosons. The quantum 
Hall phenomenology we have discussed is thus analogous to the QHE of 'spinless' (fully polarized) electrons, i.e. the 
case where the effective Zeeman gap is sufhently large that the spin degree of freedom of the electrons is frozen out. 
Let us end with a brief discussion of the more general case where internal degrees of freedom, such as spin, play a 
role. 

Polarization effects have been studied in the quantum Hall literature, and in particularly it has been shown that un- 
der certain conditions, the lowest-energy charged excitation at = 1 is a spin-textured object, a so-called skyrmion 9o| . 
Moreover, there have been studies of spin polarization effects at the edge of a quantum Hall system, predicting the 
existence of spin-textured states for sufficiently smooth confining potentials 
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92| . In the context of atomic Bose 



condensates, there are several ways of creating systems with internal degrees of freedom, promising an even richer 
phenomenology than in the quantum Hall system. One interesting approach to producing multicomponent Bose 
condensates is the simultaneous condensation of mixtures of different atomic isotopes such as *^Rb and *^Rb 
Moreover, it is possible to create spinor condensates by trapping higher-spin atoms such as ^''Rb 



95| or 23Na 
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94|. 



97| 



optical traps [98j. The advantage of this technique is that optical traps confine the atoms independently of their spin 
orientations - as opposed to traditional magnetic traps, which typically confine only one spin projection, effectively 
giving scalar condensates 



96|. 



l3 



Inspired by these experimental advances, Reijnders et al[10| have performed theoretical studies of rotating spin 
1 condensates in the lowest Landau level approximation, predicting a rich phase diagram and a number of exotic 
states. In particular, in the quantum Hall regime, they predict several series of novel non-Abelian states which are 
generalizations of the Read-Rezayi states discussed in the previous subsection. One might expect that future studies 
will continue to reveal interesting new physics in high-rotation states of multicomponent Bose condensates. 



IX. CONCLUDING REMARKS 



A summary of the research on quantum Hall physics in rotating atomic gases is necessarily preliminary, as the field 
is still highly active. At the present stage it is probably fair to say that the theoretical side is well explored - there 
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are many direct analogies to the conventional quantum Hall effect, but also physical differences, such as the expected 
occurrence of non-Abelian states in the lowest Landau level. The desirable next step would be for experiments to 
'catch up' and reach the quantum Hall regime. There is reason to be optimistic: The experimental development has 
been rapid since the first theoretical prediction of a quantum Hall effect in rotating BEC and the first experimental 

creation of a quantized vortex in the late nineties. The race towards the first experimental realization of the bosonic 
quantum Hall regime is going on at the time of writing, hand in hand with new theoretical proposals how to best 
design such an experiment. 

One can only speculate about future developments. Given the large degree of controllability of various parameters 
in cold atom experiments, one may dream of the possibility to create and manipulate anyonic quasiparticles and 
directly measure their fractional statistics, Abelian or non-Abelian. This, in turn, might be the first step towards 
implementing a topological quantum computer. But this certainly cannot be expected to happen in the near future. 
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